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Few-Body Lattice Calculations 

Martin J. Savage 

Department of Physics, University of Washington, Seattle, WA 98195-1560. 



I will discuss the recent progress toward computing few-body observables using 
numerical lattice techniques. The focus is overwhelmingly on the latest results 
from lattice QCD calculations. I present preliminary results from a lattice 
calculation of the central and tensor potentials between B-mesons in the heavy- 
quark limit. 
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i 1. Introduction 

Perhaps the greatest challenge facing those of us working in the area of 
strong interaction physics is to be able to rigorously compute the properties 
O ' and interactions of nuclei. The many decades of theoretical and experimen- 

2 ■ tal investigations in nuclear physics have, in many instances, provided a very 

precise phenomenology of the strong interactions in the non-perturbative 
regime. However, at this point in time we have little understanding of much 
of this phenomenology in terms of the underlying theory of the strong in- 



\ teractions. Quantum Chromo Dynamics (QCD). 

The ultimate goal is to be able to rigorously compute the properties 
and interactions of nuclei from QCD. This includes determining how the 
structure of nuclei depend upon the fundamental constants of nature. Any 
nuclear observable is essentially a function of only five constants, the length 
scale of the strong interactions, Aqcd , the quark masses, to„, ma, TOs, and 
the electromagnetic coupling, ae (at low energies the dependence upon the 
top, bottom and charm quarks masses is encapsulated in Aqcd) - Perhaps as 
important, we would then be in the position to reliably compute quantities 
that cannot be accessed, cither directly or indirectly, by experiment. 

The only way to rigorously compute strong-interaction quantities in the 
nonperturbative regime is with lattice QCD. One starts with the QCD La- 
grange density and performs a Monte-Carlo evaluation of Euclidean space 
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Green functions directly from the path integral. To perform such an eval- 
uation, space-time is latticized and computations are performed in a finite 
volume, at finite lattice spacing, and at this point in time, with quark 
masses that are larger than the physical quark masses. To compute any 
given quantity, contractions are performed in which the valence quarks that 
propagate on any given gauge-field configuration are "tied together". For 
simple processes such as nucleon-nucleon scattering, such contractions do 
not require significant computer time compared with lattice or propagator 
generation. However, as one explores processes involving more hadrons, the 
number of contractions grows rapidly (for a nucleus with atomic number 
A and charge Z, the number of contractions is {A + Z)\{2A — Z)l), and a 
direct lattice QCD calculation of the properties of a large nucleus is quite 
impractical simply due to the computational time required. The way to 
proceed is to establish a small number of effective theories, each of which 
have well-defined expansion parameters and can be shown to be the most 
general form consistent with the symmetries of QCD. Each theory must 
provide a complete description of nuclei over some range of atomic number. 
Calculations in two "adjacent" theories are performed for a range of atomic 
numbers for which both theories converge. One then matches coefficients 
in one EFT to the calculations in the other EFT or to the lattice, and 
thereby one can make an indirect, but rigorous connection between QCD 
and nuclei. It appears that four different matchings are required: 

(1) Lattice QCD. Lattice QCD calculations of the properties of the very 
lightest nuclei will be possible at some point in the not so distant fu- 
ture.^ Calculations for A < 4 as a function of the light-quark masses, 
would uniquely define the interactions between nucleons up to and in- 
cluding the four-body operators. Depending on the desired precision, 
one could possibly imagine calculations up to A 8. 

The chiral potentials and interactions have been determined out to 
the order where four-body interactions contribute.^ As with the EFT 
constructions in the meson-sector and single-nucleon sector, the number 
of counterterms proliferates with increasing order in the expansion, and 
at some order one looses rigorous predictive power without external 
input. Lattice QCD will make model independent determinations of 
these counterterms. 

(2) Exact Many-Body Methods. During the past decade one has seen 
remarkable progress in the calculation of nuclear properties using Green 
Function Monte-Carlo (GFMC) with the Ayis-potential (e.g. Ref.^) 
and also the No-Core Shell Model (NCSM) (e.g. Ref."*) using chiral 
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potentials. Starting with the chiral potentials, which are the most gen- 
eral interactions between nucleons consistent with QCD, one would 
calculate the properties of nuclei as a function of all the parameters in 
the chiral potentials out to some given order in the chiral expansion. 
A comparison between such calculations and lattice QCD calculations 
will determine these parameters to some level of precision. These pa- 
rameters can then be used in the calculation of nuclear properties up to 
atomic numbers A ^ 20 — 30. The computer time for these many-body 
theories suffers from the same ^ {Al)"^ blow-up that lattice QCD does, 
and for a sufficiently large nucleus, such calculations become impracti- 
cal. 

Another recent development that shows exceptional promise is the lat- 
ticization of the chiral effective field theories. This should provide 
a model-independent calculation of nuclear processes once matched to 
lattice QCD calculations. 

(3) Coupled Cluster Calculations. In order to move to larger nuclei, 
j4< 100 a technique that has shown promise is to implement a coupled- 
clusters expansion (e.g. Ref.^°). One uses the same chiral potential that 
will have been matched to lattice QCD calculations, and then performs 
a diagonalization of the nuclear Hamiltonian, after truncating the clus- 
ter expansion, which itself contains arbitrary coeflacients. The results 
of these calculations will be matched to those of the NCSM or GFMC 
for j4 ~ 20 — 30 to determine the arbitrary coefficients. This method is 
unlikely to be practical for very large atomic numbers. 

(4) Density Functional Theory (??) and Very Large Nuclei To com- 
plete the periodic table one needs to have an effective theory that 
is valid for very large nuclei and nuclear matter. A candidate that 
has received recent attention is Density Function Theory (DFT) (e.g. 
Refs.^-'^'-'^^). It remains to be seen if this is in fact a viable candidate. 
There is reason to hope that this will be useful because there is clearly 
a density expansion in large nuclei with a power-counting that is con- 
sistent with the Naive Dimensional Analysis (NDA) of Georgi and 
Manohar.-'^'' The application of DFT to large nuclei is presently the 
least rigorously developed component of this program. 

The latticized chiral theory mentioned previously can also be applied 
to the infinite nuclear matter problem. This work is still in the very 
earliest stages of exploration, but this looks promising.'^ 
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2. Lattice QCD Calculations of Single Nucleon Properties 

While the lattice QCD calculations have historically focused on the meson 
sector, both the properties of single mesons and the interactions between 
two mesons, primarily due to limitations in computational power, the last 
few years has seen an increasing number of precise calculations of the prop- 
erties of nucleons at the available pion masses. I wish to show two that will 
be of interest to the participants of this meeting. 

2.1. The Matrix Element of the Axial Current 
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Fig. 1. gA as a function of m^. This figure is reproduced with the permission of LHPC. 

The matrix element of the axial current in the nucleon, qa, is a fun- 
damental quantity in nuclear physics as it is related to the strength of the 
long-range part of the nucleon-nucleon interaction via PCAC. During the 
last year LHPC^'' has for the first time calculated qa at small enough pion 
masses where Heavy Baryon Chiral Perturbation Theory (HBxPT) should 
converge. The results of this calculation, previous calculations and the phys- 
ical value, are shown in Fig. 1. Also shown is the chiral extrapolation along 
with is uncertainty (the shaded region). The physical value lies within the 
range predicted by chiral extrapolation of the lattice calculations. 

2.2. The Neutron- Proton Mass Difference 

During the last year it was realized that one could use isospin-symmetric 
lattices to compute isospin-breaking quantities to NLO in the chiral ex- 
pansion."'^^ This is achieved by performing partially-quenched (unphysical) 
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calculations of the nucleon mass"'^^ in which the valence quark masses differ 
from the sea-quark masses (those of the configurations) and determining 
counterterms in the partially-quenched chiral Lagrangian. The NPLQCD 
collaboration found that, in the absence of electromagnetic interactions, the 
neutron-proton mass difference at the physical value of the quark masses is 

Mn-Mpl""""" =-h2.26± 0.57 ±0.42 ±0.10 MeV , (1) 
which is to be compared with estimates derived from the Cottingham sum- 
rule of Mn — -^p|physicai ~ ±2.05 ± 0.3. We see that the lattice determina- 
tion is consistent with what is found in nature. The lattice calculation is 
expected to become significantly more precise during the next year. 

3. Hadron Scattering from Lattice QCD 

To circumvent the Maiani- Testa theorem,^'' which states that one cannot 
compute Green functions at infinite volume on the lattice and recover S- 
matrix elements except at kinematic thresholds, one computes the energy- 
eigenstates of the two particle system at finite volume to extract the scatter- 
ing amplitude. The scattering amplitude of two pions in the I = 2-channel 
has been the process of choice to explore this technique, and to determine 
the reliability and systematics of lattice calculations. 

3.1. TTTT Scattering 

NPLQCD has calculated tttt scattering in the 1 = 2 channel at relatively 
low pion masses using the mixed-action technique of LHPC. Domain-wall 
valence propagators are calculated on MILC configurations containing stag- 
gered sea-quarks. The lattice calculations were performed with the Chroma 
software suite^®'^° on the high-performance computing systems at the Jef- 
ferson Laboratory (JLab). 

This calculation demonstrates the predictive capabilities of lattice QCD 
combined with low-energy EFT's. By writing the expansion of m,r<^2 as a 
function of m^//^, it has been shown that when inserting the values of 
and as calculated on the lattice. The lattice spacing effects in this mixed 
action calculation, which naively appear at are further suppressed, 

appearing only at higher orders. ^^'^^ The first non-trivial contributions from 
the partial-quenching in this calculation are shown to be numerically very 
small, and therefore a clean extraction of the counterterm in the xPT La- 
grangian is possible. The chiral extrapolation, as shown in Fig. 2, is found 
to have a smaller uncertainty at the physical point than the experimental 
value. 
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Fig. 2. The tttt scattering length in the 1 = 2 channel as a function of m'^/f^. The 
dashed curve corresponds to the unique prediction of tree-level xPT, while the shaded 
region is the fit to the results of the NPLQCD calculations.'^'^ The NPLQCD calculations 
are performed at a single lattice spacing of 6 ~ 0.125 fm. 

3.2. Ktt Scattering 

Studying the low-energy interactions between kaons and pions with K^n^ 
bound-states allows for an explicit exploration of the three-flavor structure 
of low-energy hadronic interactions, an aspect that is not directly probed 
in TTTT scattering. Experiments have been proposed by the DIRAC collab- 
oration^'' to study Kn atoms at CERN, J-PARC and GSI, the results of 
which would provide direct measurements or constraints on combinations 
of the scattering lengths. In the isospin limit, there are two isospin channels 
available to the Ktt system, J = ^ or / = |. The width of a K^tt~ atom 
depends upon the difference between scattering lengths in the two channels, 
r ~ (ai/2 ~ '13/2)^1 (where ai/2 and 03/2 are the I = and I = ^ scattering 
lengths, respectively) while the shift of the ground-state depends upon a 
different combination, AEq ^ 2ai/2 + 03/2- 

The NPLQCD collaboration calculated the K^tt^ scattering length, 
ax+TT+ 1 in the same way as we computed tttt scattering, requiring only the 
additional generation of strange quark valence propagators.^^ As af^+^+ 
was calculated (again at a single lattice spacing of 6 ~ 0.125 fm) at three 
different pion masses, a chiral extrapolation was performed. This extrapo- 
lation depends upon two counterterms, one from the crossing-even and one 
from the crossing-odd amplitudes, and the important point is that their 
coefficients have different dependence upon the meson masses. Therefore, 
both can be determined from the results of the lattice calculation, as shown 
in Fig. 3, and therefore, the scattering lengths in both isospin channels can 
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Fig. 3. The left panel shows the 95% confidence ellipses for the counterterms L5 and 
-Lxtt that contribute to Ktt scattering at NLO in xPT as extracted from the mixed- 
action lattice QCD calculation of NPLQCD at a lattice spacing of 6 ~ 0.125 fm.^^ 
The right panel shows the 95% confidence ellipses for the Ktt scattering lengths from a 
combination of a lattice QCD calculation and xPT Also shown are the 38% confidence 
ellipses from a Roy-Steiner analysis and from a 0{p'^) xPT analysis. 

be predicted, as shown in Fig. 3. This is a another demonstration of the 
combined power of lattice QCD and xPT. 

3.3. Nucleon-Nucleon Scattering 

A few years ago we reahzed^®'^^ that even though the scattering lengths 
in the nucleon-nucleon system are unnaturally large, and much larger than 
the spatial dimensions of currently available lattices, rigorous calculations 
in the NN-sector could be performed today. There are two aspects to this. 
First, it is unlikely that the scattering lengths in the NN sector are unnat- 
urally large when computed on lattices with the lightest pion masses that 
are presently available, r7i^°^> 250 MeV. Second, it is not the scattering 
length that dictates the lattice volumes that can be used in a Luscher-type 
analysis, but it is the range of the interaction, which is set by for the 
NN interaction. While the Luscher asymptotic formula are not applicable 
when the scattering length becomes comparable to the spatial dimensions, 
the complete relation is still applicable. 



p cot 6 (p) = 




- 47rA . (2) 



where S{p) is the elastic-scattering phase shift, and the sum in eq. (2) is 
over all triplets of integers j such that |j| < A and the limit A ^ 00 is 
implicit. 
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The same set of configurations that was used to produce the scattering 
lengths for tttt and Kn scattering, discussed above, were used to construct 
the NN correlation functions, and hence to extract the energy shifts between 
two nucleons at finite volume and twice the single nucleon mass at finite 
volume. At the pion masses accessible, the NN scattering lengths are found 




[MeV] [MeV] 



Fig. 4. Scattering lengths in the ISq and Si — ^-Di NN channels a sa function of pion 
mass. The experimental value of the scattering length and NDA have been used to 
constrain the extrapolation in both BBSvK^^"^" and W'^^"'''^ power-countings at NLO. 

to be of natural size, set by the inverse pion mass. Only one of the pion 
masses is within the region described by the low-energy EFT's, and as such 
a prediction of the scattering length from lattice QCD alone cannot be 
made. However, when combined with the physical values of the scattering 
lengths, an allowed region for the scattering lengths as a function of pion 
mass can be made, as shown in Fig. 4. 

4. Hadronic Potentials 

Lattice QCD cannot, in general, isolate the potential between hadrons, as 
the potential by itself is not an observable quantity. An exception to this 
statement is when both hadrons are infinite massive, and consequently their 
spatial separation can be fixed (this of course generalizes to N-hadrons). 
Currently, Silas Beane, William Detmold, Kostas Orginos and I are per- 
forming a quenched calculation of the potential between two B-mesons in 
the heavy quark limit. In this limit the light degrees of freedom (dof) of 
each B-meson decouple from the heavy-quark dynamics, and their spin, sj, 
becomes a good quantum number. The potential between the two B-mesons 
when constructed from xPT has the same form as that for NN, but with 
different values of the counterterms that enter. Therefore, calculating the 



February 9, 2008 0:48 



WSPC - Proceedings Trim Size: 9in x 6in savage ! 



9 



potential between two B-mesons may provide qualitative insight into the 
NN-system. Our calculation is not the first attempt to extract the potential 
between heavy hadrons, but is performed with a lattice spacing approxi- 
mately half of that of the previous calculations. '^'^'^^ The quark mass was 
chosen to give m-^ ^ 420 MeV and rrip ~ 700 MeV. Fig. 5 shows our prelim- 





Fig. 5. Preliminary calculations of the potential between two B-mesons in the heavy 
quark limit. The vertical axis is the difference between the central potential at a sepa- 
ration r and that at r = 7b. The fine (red) dots correspond to the leading finite-lattice 
spacing correction that must be added to the potential obtained in the lattice calculation. 



inary results for the central potentials as a function of B-meson separation 
(the tensor potential is found to be very small, and so we have also shown 
the potentials at r = \/2 and \/5 which are not separated into tensor and 
central). One clearly see short-distance repulsion in the channel with the 
quantum numbers of s-wave NN interactions, while one finds large attrac- 
tion in the other channels. There is a correction due to the finite lattice 
spacing that is to be added to the results of the lattice calculation, and 
the leading correction is shown in Fig. 5. Adding this factor gives the con- 
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tinuum potential between two infinitely massive B-mesons in the lattice 
volume, while the extrapolation to infinite volume potential has not yet 
been performed. 
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Fig. 6. Preliminary calculations of the tensor potential between two B-mesons in the 
heavy quark limit. 



The tensor interaction between hadrons has not been directly isolated in 
previous works. In Fig. 6 we show our preliminary calculation of the tensor 
potential between two infinitely massive B-mesons. We find a potential that 
is consistent with zero within the uncertainty of the calculation, indicating 
that the tensor potential in this system is no larger than Vt< 40 MeV. 

5. The Three-Body Sector 

A discussion of the three-body sector would require multiple presentations 
all by itself. I want to say just a few words about it. In order for lattice 
QCD to make connections with the three-body sector in nuclear physics 
there must be both theoretical and numerical developments. So far essen- 
tially none of the formalism has been developed with which to analyze 
lattice results. The analogous relations to the Luscher relation in the NN 
sector must be developed, and numerically explored to determine param- 
eters that can be used in lattice calculations. On the numerical side, it is 
clear from the poor quality of the NN correlation functions obtained by 
NPLQCD, that extracting any signal from the NNN sector will be very dif- 
ficult, and likely not possible on the coarse MILC lattices. We are currently 
discussing performing the contractions for the triton with the propagators 
we currently have to quantify the situation. Performing the contractions 
rapidly becomes computationally expensive. For the proton there are just 
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2 contractions, for pp there are 48 contractfons while for the triton there 
are 2880 contractions, this is simply uldL At this stage we believe that the 
extraction scattering parameters from lattice QCD in the three-body sector 
(in the near future) will require the asymmetric Clover lattices that JLab 
have applied for computer resources to produce during the next couple of 
years. 

6. Nuclear Effective Field Theory Lattice Calculations 

The last subject I wish to mention, but again I have no time to actually 
discuss it, is the recent development of performing EFT calculations on the 
lattice. I think that this is an area that has a bright future and I refer those 
that are interested to recent papers such as Refs.^'^^ 

7. Conclusions 

It is clear that lattice QCD is an important part of the future of nuclear 
physics. Lattice QCD, when combined with effective field theory, is just 
starting to make rigorous predictions of few-body observables, and we can 
expect significantly more progress as the computational resources dedicated 
to these calculations is increased. The formal tools are in place to explore 
the two-body sector, all that is required is computational power. However, 
a coherent effort involving both numerical calculations and formal develop- 
ments is presently required to move beyond the two-body sector. 

I would like to thank my collaborators in this work, Silas Beane, Paulo 
Bedaque, Kostas Orginos, William Detmold, Tom Luu, Elisabetta Pallante 
and Assumpta Parreno. 
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